LTH-418 



0^ 



Oh! 



X 



Non-zeta knots in the renormalization of the Wess-Zumino 

model? 



P.M. Ferreira & J. A. Gracey, 

Theoretical Physics Division, 

Department of Mathematical Sciences, 



Q\ • University of Liverpool, 

Peach Street, 

^ ! Liverpool, 

Q ; L69 7ZF, 

United Kingdom. 



> 

o 

(N 



Abstract. We solve the Schwinger Dyson equations of the 0{N) symmetric Wess-Zumino model 
at 0{1/N^) at the non-trivial fixed point of the d-dimensional /3-function and deduce a critical 
^ • exponent for the wave function renormalization at this order. By developing the e-expansion 

^ of the result, which agrees with known perturbation theory, we examine the distribution of 

transcendental coefficients and show that only the Riemann (^ series arises at this order in 1/A^. 
Unlike the analogous calculation at the same order in the bosonic 0{N) cp^ -theory non-zeta 
transcendentals, associated with for example the (3,4)-torus knot, cancel. 



One of the more exciting developments in renormalization theory in recent years has been 
the reahzation that there is a relation between the Riemann zeta and non-zeta transcendental 
numbers which appear in the renormalization group functions of a variety of theories and the 
positive knots of mathematics, [Q, ^, ^. Briefly, in Q it was demonstrated that a multiloop 
Feynman diagram could be easily associated with a link diagram which, after applying basic knot 
manipulations and a skeining relation, could be reduced to an underlying positive knot listed in 
the classification tables, plus trivial unknots. It turned out that one could build up a one-to-one 
correspondence between different positive knots which would appear in a Feynman diagram and 
the irrational numbers, like the Riemann zeta series, C(^)) which appeared in the renormalization 
constant of the same graph. For example, if a diagram contained the number ^(3) in its simple 
pole with respect to the regularization then it always occurred in a Feynman diagram whose 
associated link diagram reduced to the trefoil or (2, 3)-torus knot. More generally, the presence 
of C(2n — 3), for n > 3, in a renormalization constant corresponds to a (2, 2n — 3)-torus knot in the 
link diagram of the original Feynman diagram, jl], Q]. Indeed intense investigation revealed that 
other non-zeta transcendentals which are represented by independent double and triple infinite 
series were 1-1 associated with other sets of torus knots, Q. Moreover it was found that there 
was an intimate relation between the braid word representation of the knot and the structure 
of its knot number when it is expressed as an infinite sum. This empirical evidence has been 
examined with great intensity more recently and has so far remained robust, 0]. An overview of 
the status of the area can be found in ||5| . One advantage for renormalization theory of having 
such an association is that one immediately knows from the simple link diagram of a Feynman 
diagram which independent transcendentals to expect in its calculation. This can therefore be 
adapted as a calculational aid for the very high order loop diagrams in which they arise. For 
example, performing an algebraic-numerical calculation one can project onto the appropriate 
number basis for the diagram. Indeed this method was successfully used in |^] to compute all 
the primitively divergent graphs at six and seven loop which contribute to the MS /3-function of 
0^ theory. This was the first order in this theory to involve these new non-zeta transcendentals. 

Clearly to gain some insight into the underlying knot theory requires very high order mul- 
tiloop calculations which is a technical and difficult exercise, pi, 0, 0]. However, in [P] non-zeta 
transcendentals which occur in graphs with subgraph divergences were accessed in the wave 
function renormalization of 0{N) cp'^ theory through the large A^ expansion. (The knot relation 
in graphs with subgraph divergences has also been investigated in perturbation theory in [0].) 
A result had been available in |8| in d-dimensions at 0{1/N^) for the critical exponent r] which 
is related to the wave function anomalous dimension through the critical renormalization group. 
It contained the value of a two loop integral which when expanded in powers of e, where d = 
4 — 2e, contained non-zeta transcendentals. The first of these had been known for a long time 
and had been denoted in |^] by Uq2 = J2^>m>o {—^)^~"^n~^m~'^. However, it is now associated 
in the more modern approach with the (3, 4)-torus knot, Q, and hence through the braid word 
representation of this knot Uq2 has been replaced by the independent, but more systematically 
defined, non-alternating double sum F53 = Y^^>m>o n~^m~'^. This discovery yielded the expan- 
sion of the two loop integral to much higher order, Q, and established the knot association to a 
new number of loops. It also provided a basis for making a stronger connection with algebraic 
number theory, y]. Although most of the results we have mentioned have been established 
in simple scalar theories the knot approach has proved useful in gaining insight into long term 
multiloop problems in theories with symmetries. Indeed the cancellation of ((3) in the quenched 
/3-function of QED at three loops has been explained and extended to higher order cancellations 
by knot theoretic arguments in [^] . Therefore it would seem that the presence of a symmetry in 
a theory, such as gauge symmetry, can in some way suppress the appearance of these transcen- 
dental numbers. Moreover it would be interesting to understand what effect other symmetries, 
such as four dimensional supersymmetry, will have on the structure of the renormalization group 



functions. This is the subject of this paper where we will study the Wess-Zumino model, [pH] , 
with an 0{N) symmetry to 0{1/N^) in a large N expansion. The aim is to first calculate the 
critical exponent, rj, corresponding to the wave function renormalization and then to examine the 
expression obtained for it to ascertain how super symmetry modifies the appearance of zeta and 
non-zeta transcendentals compared with the underlying bosonic (/)^-theory. This would there- 
fore extend some of the analysis of Q. We note that in the underlying bosonic model, 0{N) 
c/)^ theory, it was possible to establish not only the first location of Uq2 in the wave function 



renormalization at 0{1/N^) but also to provide its explicit coefficient, |^, 12|. Knowledge of 
this is clearly important if one is ever to have a systematic knot based method to perform the 
renormalization since the coefficients of the transcendentals that appear in some renormalization 
scheme will also need to be determined. Therefore the large N results will provide important, 
non-trivial and independent checks. We recall that large N calculations in models such as the 
Wess-Zumino model are possible at 0{1/N^) through critical renormalization group methods 
and solution of the Schwinger Dyson equations order by order in 1/N at a non-trivial fixed 
point of the d-dimensional /3-function. The method had been originally applied to the 0{N) 
bosonic a model in an impressive series of articles, [^, Hj. More recently, similar calculations 
in supersymmetric theories have been made possible in [^ by extending the methods of [|l^] 
to superspace for the 0{N) Wess-Zumino model which, it is hoped, can be extended to study 
supersymmetric gauge theories. 

We now turn to the specifics of calculating r] in the Wess-Zumino model with an 0{N) 
symmetry and note that its action in superspace is given by 

S = I d'^x Ud^e f¥<^^ + ^\ - ]^j dHa^^ - ]^j <ma^A (1) 

where $* and a are chiral superfields, 1 < i < A^ and g is the coupling constant which has 
been rescaled into the a kinetic term to ensure that the interaction in d-dimensions is in a 
form for which the conformal integration technique or uniqueness, [|l^], can be readily applied 
to ease the calculation of high order integrals. The four loop d-dimensional /?- function of the 
model is deduced from the result for a general symmetry group given in ||l5| which used the 



non-renormalization theorem, [[ll| , |l^, and generalized earlier results, |17, |18[. Specifying to an 
0{N) group it is given by, |lF 



m = \{d-A)g + ^(iV + %' - 2{N + l)g' + Q(iv2 + lliV + 4) + 3(iV + 4)C(3) ) 5' 
+ (-{N^ - 36A^2 - 84A^ - 20) - ^N"^ + 16iV + 24)C(3) 



+ -{N + 4)^C(4) - 10(iV + 2){N + 4)C(5) j 5' + 0{g^) (2) 

where C(n) is the Riemann zeta function. As we will be calculating information relevant to the 
$-superfield renormalization, we note also that, |15|, 



7*(5) = 9 - \{N + 2)g^ - ^ (iV^ - lOiV - 4 - 24C(3)) 5' 



— (3iV3 _^ iQjsi"^ _^ 152^ + 40) _ ii^(Ar2 _ 4Ar _ 36) (at + 4) 
- 3C(4)(7V + 4) + 20C(5)(iV + 2) ) / + 0(/) (3) 

The d-dimensional fixed point where we will analyse the Schwinger Dyson equation is given by 
the non-trivial zero of the /^-function, (y). Specifically expanding in powers of 1/A^ it is given 



by, il, 

+ f 32e - 176e2 - 8(6C(3) - 37)e^ 



+ y (60C(5) - 9C(4) + 18C(3) - 4)6^ + 0{e')^ ^3+0 (^) (4) 

Evaluating the renormalization group functions at gc defines the critical exponents. So, for 
example, rj = 'y^{gc)/2. In [|l^] an expression for r] was deduced at 0{1/N'^) by replacing the 
lines of the Schwinger Dyson equation, truncated at the same order, by the asymptotic scaling 
forms near Qc of the propagators of the respective fields. For an Euclidean space-time these are, 



(p2)M-" 

B6^e - e') 

(p2 



{a{-p,e)a{p,e')) ~ : j;_; ' (5) 



as p ^ CX3, where we have set d = 2fi, A and B are unknown momentum independent amplitudes 
and 

a = /i-l + i??,/3 = l-7/ (6) 

are the exponents of the <1>* and a superfields. Their canonical dimensions are fixed by ensuring 
that each term in the d-dimensional superspace action is consistent with S being dimensionless. 
Due to the non-renormalization theorem, |11, 16 1, there is no vertex renormalization and its 
associated exponent, ordinarily denoted by X: is zero to all orders in 1/A^. In particular in the 
critical region the $ equation of fig. 1 becomes 

= a{a — fi + 1) + a{a)z (7) 

where z = a?{fi — a)a{fi — I3)A^B and a{x) = T{fi — x)/T{x). The a equation is 

= a(/3-^ + l) + ^Nza{l3) (8) 

where we have chosen to express each representation of fig. 1 in coordinate space through a 
Fourier transform. Thus eliminating z gives 

^ 4a(a)a(/3 - /i + 1) ,g. 

Due to the chirality property of the superfields no new two and three loop graphs occur at 
0{1/N'^) unlike the bosonic model and therefore one merely expands out (Ph to 0{1/N'^) to 



deduce 772, where ij = X^i^i Vi/^^- Thus, |14] 

4r(2/i - 2) 

r(;u)r2(^-i)r(2-/i) 



Vi 



V2 



iP{2-f,) + ^(2^-2) - ^(^-1) - ^(1) + ^ 



vt (10) 



2(/x- 

where ip{x) is the logarithmic derivative of the Euler F- function. 

The absence of these higher order graphs and the finiteness of the vertex means that one can 
use the Schwinger Dyson approach to compute 773 which we believe is the first instance this has 
been attempted. Ordinarily one has to use the conformal bootstrap method, [|19| , |8[, to do this. 



For the Wess-Zumino model, however, it does not seem possible to apply it since that method 
relies on the vertex dimension being non-zero. Moreover the chirality of the fields ensures that 
there are only two new topologies at 0{1/N^) and these are given for the respective Schwinger 
Dyson equations in figs 2 and 3. If we denote the values of the graphs for the <1> equation by Sj 
and those for the a equation by IIj then the respective Dyson equations to 0{1/N^) are now 

= "(" ~ ^ + 1) + z + z^a^a)a^(f3)a(fi - a - /3)J:i + Nz^a\a)a\/3)a(fi-a-m2 
a[a) 

= ^""((^ - /" + '^) + z + z^a%a)a\p)a{fi-2a)Ui + Nz^a\a)a^{P)a{fi - 2a)U^ (11) 

Na{jj) 

We note that although the factor a(/i — a — /5) that appears with each of the higher order 
corrections in the ^ Dyson equation is singular when expanded in powers of 1/N it will turn out 
that in the computation of the diagrams in momentum space there will be a compensating factor 
of a{a + /3) in its final value. Therefore the final overall contribution will be free of infinities. 

The finiteness of the vertex means that one does not need to perform the usual renormaliza- 
tion of the Schwinger Dyson equations. Therefore it only remains to compute the four graphs. 
As is usual in superfield calculations the first step in such an evaluation is to carry out the ma- 
nipulation of the supercovariant derivatives present in the super- Feynman rules which is known 
as the D-algebra. (For an introduction see, for example, |^0|.) For the topologies illustrated in 
figs 2 and 3 this will result in the distribution of D-operators in momentum space on various 
lines. For the non-planar topology the two upper lines joining to the external vertex are each 
modified by a D-operator. For the four loop topology the corresponding lines also gain the same 
operator as well as the completely internal horizontal propagator. The effect of a D-operator is 
to reduce the power of the exponent of that line by unity. For a (p line this means the power of 
its propagator will vanish whereas that of a a line will be non-zero. One might expect that for 
the graphs of the ^-equation the values of the graphs where the D-operators are distributed on 
the lower propagators are different. Explicit calculation shows that, at least to the order in 1/A^ 
which we are interested in, this is not the case. Indeed it turns out that these four graphs are 
straightforward to compute by the rules of conformal integration. The two non-planar graphs 
involve the same basic two loop integral which we denote by -P(^). It is defined graphically 
in fig. 4 in the coordinate space representation where one integrates over the location of the 
vertices. Thus, 

ni = ^^^(/^) (12) 

and 

a(/i-a-/3)Si = P{fi) (13) 

where we have included the compensating factor a{a + (3) on the left side of the value of the 
graph. This integral, P^fJ-), has not appeared in previous 0(1/ N^) calculations in other models 
and if we define 

then its e-expansion, where d = 4 — 2e, is known to very high orders, M, E^]. In particular 

n(^) = 6C(3) + 9C(4)e + 7C(5)e2 + ^[({6) - 2e{3)]e^ - i[91C(7) + 120C(4)C(3)]e4 

- -^[5517C(8) - 512C(5)C(3) - 768Ue2]e' + 0{e^) (15) 

The presence of the non-zeta transcendental at O(e^) encourages us to believe there will be a 
(3,4)-torus knot in the wave function renormalization at 0(1/A^^) at some very high order in 



perturbation theory when one disentangles the information encoded in the relation t] = ^^{gc)/2 
using the 0{1/N'^) corrections to Qc available in the /3-function critical exponent 0^2 computed in 
|14|. For 112 the integral in fact occurs in the /^-function calculation of the Wess-Zumino model, 
1 14 1, as well as in the original work of [|l^]. Therefore we merely quote the final result 



n. 



a^(l)a^(2^-2) 



i>\^i) ^ !.(/.) 



+ 



+ 



(/i-2) 



3e(^) 



ril 



(16) 



where G(//) = V'(/" - 1) - V''(l), ^(/^) = ^-'(2/" - 3) - i)'{2, - /i) - V'(/" - 1) + V''(l) and 
^{fj) = il>{2iJi — 3) + ^(3 — fi) — ip{fi — 1) — ip{l). Therefore, like Si and Hi, 112 only has 
zeta transcendent als in its e-expansion. For S2 elementary manipulations yield the three loop 
integral of fig. 5 which is illustrated in the coordinate space representation in the case when 5 = 
0. It is evaluated by integration by parts using a temporary analytic regularization introduced 
by shifting the exponents of some of the lines by an infinitesimal quantity, 5. This is necessary 
because the set of integrals which result are divergent though their sum is finite. The symmetric 
choice illustrated in fig. 5 is to minimize the number of subsequent graphs that need to be 
calculated as well as to ensure that these can still be computed by conformal techniques. We 
find that 



a{n 



a 



13)^2 = a{2fi-2) 



l^i^) 



+ 2Q\^l) 



[n) 



where ^^(/i) = ip"'{iJ, — 1) — ^'"{1). Hence we can now substitute for the values of the integrals 
in (|Tl| ) and expand out the leading term of @ to 0{\/N^) to deduce that in d-dimensions 



m 



v! 



(7,^ -26, + 25).,^^^^ (3,^ -10, + 9).^^^ 



4(/i-2)2 
5(,-l)2 



12 



n{fi) + 



4(/x-2)2 
(26,^ - 174,3 + 435,2 



(,-l)2e2(,) 
479, + 195) 



(7,2 - 16, + 10) 
4(,-2)2 



e(,) + 



2(2,-3)(,-l)(,-2)3 
(44,4 _ 257^3 _^ 557^2 _ 53^^ _^ 
2(2,-3)2(,-l)2(,-2)2 



vI/(,) 



(18) 



In the final expression it turns out that the integral -P(,) has cancelled. This is a direct result of 
the -D-algebra, and therefore supersymmetry, which modifies the exponents of each integral to 
be similar. A cancellation of some sort between Si and Hi was to have been expected, however, 
from knowledge of the e-expansion of the expression obtained for 773 by ignoring the higher 
order contributions Sj and IIj. In other words that obtained by iterating (0) to 0{1/N^). This 
approximation, known as Hartree Fock, in fact agrees with the known four loop perturbative 
result, (pp, which implies that although Si and Hi are non-zero, their ^(3) contribution at three 
loops must at least cancel as well as C(3) and C(4) at four loops. (We note the Hartree Fock 
contribution to r/4 also agrees with the four loop perturbative result.) It could have been the 
case that Si and Hi took different values but had e-expansion agreement at low orders. Indeed 
this is the situation at four loops for S2 and n2 whose low order e-expansions are the same, 
correctly involving ("(5) as the first term but differing first at level C(7). They give consistency of 
the Hartree Fock 773 to four loops but new non-trivial contributions will enter at five and higher 
loops. Therefore the explicit cancellation of the -P(,) function is novel and indicates that unlike 
its underlying bosonic 0{N) cp'^ theory, there will be no non-zeta knots at all orders in the wave 
function renormalization at 0{1/N^). Thus the presence of four dimensional supersymmetry 
would appear to reduce the number of complicated knots which will arise. This is in contrast 
to the case of the 0{N) supersymmetric a model in two dimensions where it is expected that 
its T/3 retains the integral yielding a (3, 4)-torus knot, |2^, |6|. One can understand this better 
by considering the simple properties of supersymmetric theories in differing dimensions. For 



instance, a four dimensional M = 1 supersyninietric scalar theory involves chiral superfields but 
in two dimensions it is theories with M = 2 supersymmetry which are chiral and not TV = 1 
theories such as the 0{N) supersymmetric a model. Therefore this would suggest that theories 
with A/" = 2 in two dimensions should not have non-zeta knots at 0{1/N^). 

Although (^) clearly will agree with the 4-loop perturbative result for 7$ (g) since the higher 
order graphs give non-zero contributions first at five loops, one question arises as to whether the 
correct signs and weighting factors have been included. However, the topologies of figs 2 and 3 
arise in the computation of the critical exponent which determines the /3- function at 0(1 /N'^), 



|14|. Therefore we have been careful to include the corrections to the Schwinger Dyson equations 



with the same relevant factors as in |14], since these were essential to gain agreement with the 
4-loop perturbative /J- function, (|2|). In light of these remarks we can now deduce several of the 
higher order coefficients of 7$(5). If we denote the 0{1/N^) part of "f^{g) as 

00 
7*(<?) = g + Y.^CrN^ + drN + er)N'"''g'- (19) 

with 62 = then it is easy to deduce that 
C5 = ^[3C(4)-2C(3)-1] 

4 = ^[18C(4)-30C(3)-29] (20) 

1^ ^- ^^^.2. 



65 = — [-450C(6) + 1416C(5)-171C(4) + 102C(3)-108C(3)+304] 



and 



C6 = ^[6C(5)-3C(4)-2C(3)-1] 

4 = ^[80C(5)-162C(4) + 248C(3)-81] (21) 

ee = ^[ - 230C(7) + 1140C(6) - 906C(5) - 36C(4) - 5C(3) - 108C(3)C(4) + 414(^^(3) + 14] 

Therefore there are only two unknown coefficients to be determined in the polynomial of N to 
complete the five loop result. Interestingly the first contribution to the perturbative coefficients 
from the higher order four loop diagrams of figs 2 and 3 is the C,{7) term at six loops. Although 
we have expressed our result in d-dimensions we cannot quote a value for 7/3 in three dimensions 
since the expression diverges there. This is on a par with the critical exponent u} = — P'{gc) 
which is also singular in three dimensions but at 0{1/N'^), []T4| . 

In conclusion although we have ruled out non-zeta knots at 0(1/A^'^) in the wave function 
renormalization of the 0{N) Wess-Zumino model it is not clear if this would be the case at all 
orders in perturbation theory or for the other renormalization group functions. For instance, 
one way of ascertaining whether the (3, 4)-torus knot survives at six loops in, say, the /3-function 
would be to analyse those Feynman diagrams with no subgraph divergences which contributed 
in (/>^ theory to its /?- function, Q, and determine their value in the context of the Wess-Zumino 
model in superspace. Indeed it may be the case that all of these topologies are excluded at the 
first stage by the chirality of the fields. Alternatively if that is not the case then it would be 
an intriguing exercise to see how the -D-algebra modifies the sum of the contributions from each 
diagram and if there is an overall cancellation. 
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Fig. 1. Leading order Schwinger Dyson equations. 





Fig. 2. Additional graphs contributing to the $ Dyson equation at 0{1/N^). 





HAAAA/ 



Fig. 3. Additional graphs contributing to the a Dyson equation at 0{1/N^). 




Fig. 4. Definition of the integral P{^). 




Fig. 5. Intermediate integral in the calculation of Sj 



